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ABSTRACT 

Precision weak gravitational lensing experiments require measurements of galaxy 
shapes accurate to < 1 part in 1000. We investigate measurement biases, noted by 
Voigt and Bridle (2009) and Melchior et al. (2009), that are common to shape mea- 
surement methodologies that rely upon fitting elliptical-isophote galaxy models to 
observed data. The first bias arises when the true galaxy shapes do not match the 
models being fit. We show that this "underfitting bias" is due, at root, to these meth- 
ods' attempts to use information at high spatial frequencies that has been destroyed 
by the convolution with the point-spread function (PSF) and/or by sampling. We 
propose a new shape-measurement technique that is explicitly confined to observable 
regions of fc-space. A second bias arises for galaxies whose ellipticity varies with ra- 
dius. For most shape-measurement methods, such galaxies are subject to "ellipticity 
gradient bias." We show how to reduce such biases by factors of 20-100 within the new 
shape-measurement method. The resulting shear estimator has multiplicative errors 
< 1 part in 10^ for high-S/TV images, even for highly asymmetric galaxies. Without 
any training or recalibration, the new method obtains Q = 3000 in the GREAT08 
Challenge of blind shear reconstruction on low-noise galaxies, several times better 
than any previous method. 
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The appearance of background galaxies are subtly distorted 
due to weak gravitational lensing (WL) by mass along the 
line of sight. This WL effect is detectable as a breaking of 
symmetries in the intrinsic background scene, in particular 
a shearing that breaks the isotropy of the galaxy ellipticity 
distribution. Because this "cosmic shear" is a direct mea- 
sure of the metric fluctuations in the Universe, it is readily 
used to constrain the dark matter and energy that dominate 
the gravitational behavior of our U niverse. A review of the 
cosmic shear method is provided bv lHoekstra fc JainI (|2008l ) 
and the role of WL in the con t ext of dark energy studies is 
summarized in lAlbrecht et all (|2006l ). 

The cosmic shear signal is subtle and difficult to mea- 
sure to high accuracy. The shear along a typical line of sight 
on cosmological scales causes a ~ 2% change in axis ra- 
tio of a background source. Propagation of shear-estimation 
errors into cosmological parameter inferences suggests that 
systematic errors in shear estimation must be held below 
~ 1 part in 1000 of the expected shear to avoid significant 
degradation of cosmological inferences (|Huterer et aZ.|[2006l : 



lAmara fc RefregieJ |2008| ). Roughly speaking this means 
that systematic errors in the ellipticity e of individual galax- 
ies must be held below 1 part in 1000, i.e. {Se)/e < 10~^. In 
comm unity-wide blind tests of shear recovery from simulated 
data iHevmans et al.ll2006l : iBridle et al.ll2009bl '). no meth- 
ods have yet achieved this level of accuracy, even for sim- 
ulated simplified images with unrealistically high signal-to- 
noise (S/N) levels and perfect knowledge of the point spread 
functio n (PSF). As most recently reported bv lBridle et al.l 
(|2009bl ). methods that could be practically applied to real 
survey data produced shear estimates that were in error by 
~ 1% or more of the typical cosmic shear signal. 

A biased and/or noisy shape measurement method can 
nonetheless be applied to real data if its biases are calibrated 
and compensated. A simple scheme is to fit the measured 
shear 7meas to the true (simulated) shear 7truG via 



= (1 -I- m)7truc + c. 



(1) 



The multiplicative and additive error terms m and c derived 
from simulations could then be applied to measured dataQ 
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^ We write equations treating 7 as a scalar, for simplicity. In 
practice the shear, ellipticities, additive and multiplicative error 
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The difficulty is that m and c will depend upon sets of pa- 
rameters TTg and TT* describing the properties of the source 
galaxies and point-spread function (PSF), respectively. If 
the simulation does not properly mimic the real popula- 
tions, then the empirical correction will be in error — and we 
may not even know which parameters are the critical ones 
to emulate properly. A Bayesian formalism can be used to 
derive an unbiased shear estimator 7 from a collection of po- 
tentially biased shape estimates if we know all three of: 
the intrinsic distribution of galaxy shapes P(e); the change 
dP{e)/d'y of the shape distribution under applied shear; and 
the c onditional probab ility P{e\e) of the measurement pro- 
cess (|Miller et al\ \20m). T he need for s u ch kn owledge was 
previously noted bv .Bernstein fc JarvisI (|2002l . BJ02). The 
calibration is difficult, however, because the measurement 
process noise and bias will depend upon galaxy and PSF 
properties, P(e|e, TTg, tt*), so we need to know this function 
and also the distribution P{e,Hg) of all galaxy subsets of 
interest. Characterizing all of these functions by simulation 
will be hard because of the wide variety of galaxy shapes 
on the sky and PSF variation, especially if we do not know 
exactly what parameters are important. 

It is therefore of great interest to produce a shape- 
measurement algorithm that is minimally biased, and which 
produces distributions P{e\e) that are robust, i.e. charac- 
terized by a few well understood properties of the galaxy or 
PSF. Any corrections that must be derived empirically (via 
simulation) should be as small as possible to minimize errors 
that arise from differences between the simulations and the 
real observations. 

Furthermore the calibration depends on knowing the re- 
sponse dP{e)/d^ of the noise-free, unsheared ellipticity pop- 
ulation P(e) to an applied shear. If the shape-measurement 
methodology misestimates dP/dy even in the absence of 
noise, then the calibration will be in error even if P{e\e) 
is known. If dP/d'y depends upon the detailed parameters 
■Kg and/or tt*, then the calibration will depend upon knowl- 
edge of the complete distribution P{e,ng,n-i,) at all points 
in the survey. 

In this context it is worrisome that IVoigt fc Bridi3 
l|2009l, VB) and lMelchior et al\ (|2009l . MBLB) demonstrate 
that shape-measurement techniques based on model-fitting 
to the observed data can be biased by > 1 part in 100 even in 
the absence of measurement noise. These biases arise when 
the models being fit are not well matched to the galaxies 
being measured, for example if an exponential-disk model 
or Gauss-Laguerre expansion is being fit to a deVaucouleurs 
profile galaxy. In principle the shape measurement bias can 
be reduced by adding degrees of freedom to the model so 
that it spans all possible galaxy shapes, for example by using 
a complete set of basis functions to model the galaxy. This 
will, however, come at the expense of increased uncertainty 
on the inferred shear as the model fits become increasingly 
degenerate. 

In this paper we examine these "underfitting" biases 
more closely, and devise a shape-measurement technique 
that avoids the issue. Then we investigate a more general 
effect revealed by the VB tests, in which a galaxy that has 



coefficients all have two components. We will also ignore the dis- 
tinction between the shear 7 and the reduced shear g = 7/(1 — k). 



ellipticity varying with radius (or more precisely, with the 
scale of the weight function) leads to a mis-estimation of 
the applied lensing shear. We will show how this "elliptic- 
ity gradient bias" can be quantified and corrected using our 
newly proposed shear- measurement method. After control- 
ling these two biases we are able to measure shear from 
high-S/N galaxies to the desired accuracy. 



2 ROUNDNESS-TEST METHODS 

BJ02 introduced a geometric definition of galaxy shape 
which is intended to eliminate the need for empirical shear 
calibration factors, at least for high-S/TV objects. This shape 
definition renders the derivative dP{e)/dj a universal func- 
tion of P{e), independent of any galaxy or PSF parameters. 
We briefiy summarize the notion of "geometric" shape as- 
signment. 

The shearing action of lensing that we seek to detect 
is a mapping of the unlensed galaxy appearance lu to the 
apparent surface brightness /a(x) via /a(x) = /t,(S^^x). We 
take a unit-determinant version of the shear, so that for a 
shear along the a:: axis, 

For shear at an angle /3 to the x axis, we apply the ap- 
propriate rotation matrices to S and define r^i = cos 2/3, 
?72 ~ ri sin 2/3. It is more common to write a traceless mag- 
nification matrix 

The two variables for shear are related by 7 = tanhr;/2. 
A shear-measurement method must infer the components of 
the applied shear rj from a collection a of shapes assigned 
to galaxies i £ 1, . . . ,N. One crucial question is: how do we 
define for a galaxy with some /^(x)? 

We will call a galaxy "perfectly elliptical" if all of its 
isophotes are true ellipses with constant axis ratio a/b and 
position angle p. For perfectly elliptical galaxies, the galaxy 
e is straightforwardly defined as the isophotal ellipticity 
(a^ — 6^)/(a^ -|-fe^). Furthermore an applied shear t] maps ev- 
ery such galaxy into another perfectly elliptical galaxy with 
ellipticity e', and the map e'{e,ri) has the desired property 
of being independent of all other galaxy parameters. For an 
applied shear rj <^ 1 and an unlensed galaxy ellipticity ori- 
ented along the x axis, the map can be linearized to (B J02) : 

e'l = e+{l-e^)r,cos2l3 

e'2 = 7? sin 2/3 ^ ' 

Hence if all galaxies were perfectly elliptical and it were 
possible to observe their pre-seeing shapes /a(x) at infinite 
S/N, then a perfect shear estimator would arise from simply 
measuring the ellipticity of any single isophote. 

Among the many problems with the above method 
is that galaxies do not have perfectly similar elliptical 
isophotes. BJ02 show that one can assign an ellipticity e to 
a generic galaxy /a(x) such that it is transformed by shear 
under the same model-independent formulae as the perfect 
ellipticities above. To do so we define a pair of roundness 
tests ti{Ia),t2{Ia) for which: 
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(i) There is one and only one shear rjt for which 

tl[/a(S^,x)] =f2[/a(S^,x)] =0. (5) 

(ii) If ti{Ia) ~ t-ii^Ia) = holds for some la, then it also 
holds after any rotation Rg is applied to la- 

(iii) The form of the roundness test is invariant under 
shear of the target galaxy. 

The first two conditions mean that a galaxy which nulls 
ti and t2 can be considered "round" in the sense that any 
rotated version of the galaxy is still round, and any sheared 
version of the galaxy is no longer round. If a chosen galaxy 
nulls ti and t2 at test shear T^t, we assign the galaxy the 
ellipticity e = tanh(r;t), since it is a "round" galaxy sheared 
by r)t- The third condition then guarantees that this assigned 
ellipticity will transform under shear as in Equation Q. 

It is often desirable for the form of ti to depend upon 
properties iig and tt* of the galaxy and PSF, e.g. adjusted 
for galaxy and PSF sizes to improve the measurement S/N 
level. We will emphasize later however that it is critical that 
the roundness test be the same for sheared and unsheared 
versions of the same galajcy. We will assume that the round- 
ness test is a linear function of /a(x), and hence is also a 
linear function of the Fourier-domain intensity /a(k). It is 
typical to adopt a roundness test that has a quadrupole 
form, since shear induces a quadrupole on a circular galaxy. 
A radial weight function w remains free: 



tl{Ia) 



t2{Ia) 



dxdy Ia{S^x)w{x + y ){2xy) 



(6) 



The result of the roundness tests will usually depend 
upon the choice of origin for the coordinate system — 
although there are methods for which this is not true, which 
work on the autocorrelation function or power spectrum 
of la llvan Waerbeke et all Il997l : iHosseini fc Bethgel l2009l : 
[Zh andboosh . To fix the coordinate origin xo we can cre- 



ate two additional null tests. These tests must satisfy the 
three conditions given for the shear null test, with the dif- 
ference that they are null for a unique translation xo rather 
than a unique shear. The obvious choices for linear centroid 
null tests are dipole moments with some radial weight Wx'- 



tx{Ia) 



ty{Ia) 



dxdy Ia[Sr,{x - Xo)]Wa:{x^ + y^)x = 0, (7) 



dxdy Ia[S,,{x- ■x.o)]wx{x + y )y = 0. 



The centroid and ellipticity are hence determined by search- 
ing the {xojyojTiitV^} space for the point where tx,ty,t\, 
and t2 are all nulled. When all null tests are linear functions 
of la, the process can be equivalently cast in either real or 
Fourier space. 



2.1 PSF correction with basis-set methods 

The observed galaxy image /o(x) = /a(x) (g) T(x) is always 
the lensed galaxy image convolved with some instrumental 
transfer function (PSF) T. We must however access the pre- 
seeing function la somehow to apply the roundness tests. 
The technique proposed by BJ02 is to assume that the 



roundness (and centroid) tests are all weighted integrals over 
la with weight functions that are members of a set V'i(^) of 
functions that form a complete basis for functions on the 
plane. If the axe orthonormal and complete, then there is 
a unique set of coefficients hi such that 



/a(S„x) = 



^(x) 



d X IaiSr^X.)lpt{x) 



[^/'»(S-'x)®T(x) 



(8) 



(9) 
(10) 



The assignment of galaxy shape e thus proceeds as follows: 

(i) For some trial rj (and centroid choice), shift and shear 
the basis functions ipi and convolve them with the PSF. 

(ii) Using the last of the equations above, execute a linear 
least-squares fit to the observed galaxy to determine the 
coefficients hi. In the absence of noise, the hi that provide 
the decomposition in the first line must be a solution of this 
fit with ~ 0, hence a best fit. 

(iii) Iterate the trial shear and centroid to null the four bi 
which embody the roundness and centroid tests. From the 
second equation we see that these bi would be the result of 
applying the centroid test to the pre-seeing image la- The 
galaxy is assigned an ellipticity corresponding to the shear 
that nulls the tests. [More generally the null tests are four 
functions of the bi, but the principle remains similar.] 

An obvious drawback to this method is that the basis 
set must be infinite to span the full variety of all the galaxies 
on the sky. Real data will not yield a unique solution to an 
infinite coefficient vector {bi} because the data are sampled 
on a bounded region. Furthermore there must be regions 
where the transform r(k) of the PSF is identically zero, 
which will induce degeneracies into the fit of an infinite basis 
set, even when there is no noise in the observations. We must 
ther efore apply t he method with a finite, incomplete basts 
set. iLewid (|2009l ) notes that fitting methods can become 
biased estimators of shear when the models do not fit the 
galaxies exactly. Here we explore the conditions under which 
biases develop. 



2.2 Bias from incomplete bases 

We simplify the situation by assuming that there is one co- 
efficient R — bi that serves as a roundness test, and another 
coefficient H = bj that we will consider "hidden" because we 
will later truncate it from the model. Assume that the galaxy 
is round in some coordinate system, which means that a 
complete decomposition of la would yield values iitruc = 
and -fftruc. As illustrated in Figure [1] the the fitting process 
may induce non-zero correlation between R and H. 

Next consider that truncating the fit is equivalent to 
setting _ff = and seeking the best fit for R. If Cov(i?, H) 7^ 
and //true 7^ 0, then the best-fit value R will be non-zero. 
This will be manifested as an error m the assigned ellipticity. 
Note that the uncertainty on R will be larger when H is free 
than when H is truncated; in cases with finite noise, this may 
eventually produce some upper limit to the number of basis 
functions that can be included in the model. Hence when 
the measurement process correlates the roundness test with 
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Figure 1. Illustration of bias induced by underfitting: in this sim- 
plified picture, we plot the contours of for a fit to a galaxy with 
a model having two degrees of freedom, R and H. R is a, "round- 
ness test", i.e. an indicator of shape, while H is some additional 
parameter needed to obtain a complete model of the galaxy. The 
true values -Rtrue = and -fftrue 7^ are obtained from a full fit. 
However if we truncate the model so that H is fixed to zero, then 
the best fit moves to R 7^ 0, producing a bias in the measured 
shape. Fitting to the full {R, H) model may be disfavored, or even 
impossible to execute, because the resultant uncertainty on R can 
become larger. However the truncated model will be biased. 

other members of the basis set, there is an inevitable tradeoff 
between accuracy and bias in this fitting method. 

In an ideal case of basis- function fitting, (1) the noise in 
the image is white, (2) there is no PSF, and (3) the data is 
available over the full plane. In this case the coefficients of 
an orthonormal basis will be uncorrelated, so Cov{H, R) = 
and there will be no bias from fitting a truncated basis set. 
Correlations between basis coefficients are induced, however, 
by finite sampling, bounded data range, and especia lly by 
convolution with the PSF (BJ02, iBerrv et al.] 120041 '). The 
roundness test is typically a quadrupole, with m = ±2 az- 
imuthal symmetry, in the sheared coordinate frame where 
the object is round. Introduction of a covariance between 
R and a hidden parameter with azimuthal order tuh re- 
quires that the PSF have elements with azimuthal order 
m,t = ±mH ± 2. "Selection rules" for covariances induced 
by finite sampling and data range are more complex. 

We can tolerate biases in the measurements of individ- 
ual galaxy ellipticities if we can be assured that they will go 
to zero when averaged over the full (isotropic) population of 
pre-lensing galaxy shapes. For the shear derived from a full 
population to be biased due to truncation of some parameter 
H, we need 

{HtrucCoY{H,R))^0. (11) 

Both //true and the covariance must be evaluated in the 
sheared coordinate system in which the galaxy is round. 
Consider the case when H is the coefficient of a 



monopole basis function with high radial order, which is 
truncated from the fit to the data. For example VB attempt 
to fit a deVaucouleurs-profile elliptical galaxy with a basis 
set consisting of one (or a few) Gaussian ellipsoids. In the 
frame where these galax;ies appear round, the galaxies have 
a monopole residual to the Gaussian fit. If the PSF is circu- 
lar in the observed frame, then it is elliptical in this sheared 
frame, with an ellipticity 90° out of phase with the original 
galaxy ellipticity. The m = 2 components of this (sheared) 
PSF will induce correlations between the amplitude H of 
the truncated monopole basis function and the quadrupole 
component of the fit. These induced biases will always be 
aligned with the galaxy ellipticity. The result will be a mul- 
tiplicative shear calibration error. 

Underfitting can also result from unmodelled galaxy 
terms beyond the monopole. For example if galaxies tend to 
be disky, they have m = ±4 components with a well-defined 
phase relation to the galaxy ellipticity. These terms are cou- 
pled to the roundness test by the PSF, potentially causing 
systematic over- or under-estimation of the ellipticity. 

A circular PSF acting on an isotropic population of 
galaxies has no preferred direction and hence cannot induce 
an additive shear error. If, however, the PSF has anisotropic 
components, particularly at m = 2, then underfitting can in- 
duce additive errors by correlating the quadrupole roundness 
test with unmodelled monopole components of the galaxies. 



2.3 Circular- model methods 

Several shear-measurement methods assign ellipticities to 
galaxies by fitting the observed data with models of per- 
fectly elliptical galaxies. We note that these are subject to 
the same underfitting biases. These methods work as follows: 

(i) Select a circularly symmetric model galaxy function 
'0o(x). The function may have additional parameters for size, 
etc., but this is secondary to our discussion. 

(ii) For a trial centroid xo, shape rj, and fiux scaling /, 
produce the model observed image i/'(x) ~ /V'o[S^^(x — 

xo)l (g) r. 

(iii) Calculate the likelihood {e.g. x^) of the observed data 
matching the sheared model. 

(iv) Find the shear, centroid, and fiux that maximize the 
likelihood of the data. The galaxy is assigned an elliptic- 
ity equal to the model shear that maximizes the likelihood. 
Some methods choose to estimate the likelihood-weighted 
mean values rather than the maximum-likelihood values, 
but for the high-S/N cases we consider now, the difference 
is unimportant. 

iKuiikenI (|l999h and IVoigt fc Bridi3 (|2009l ) use a method in 
which the circular model is a sum of Gaussians; iMiller et al\ 
(! 2007l ) have an exponential or deVaucouleurs profile circular 
model. 

The method will clearly succeed without bias (for high 
S/N) when the galaxies are indeed sheared versions of the 
proffered models. However for more general galax;ies, these 
methods will be subject to the same biases as the truncated- 
basis-set methods described above. Consider the set of five 
functions 

Ipo, 
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^2 



9^0 ^ _ dlpo 

dx ^ ^ dy ^ 

a^o[s-i(x)] 

drii 

av.o[s-^(x)] 



(12) 



such that we can approximate 

/^o[S^^(x-xo)] ^ /oi/'o(x) 

+ (/oa;o)V'x(x) + {foyo)ipy{x) 

+ (/or?i)^i(x) + (/o772)v)2(x). (13) 

for small departures from a centered, circular model. The 
circular-model-fitting technique will consider xo = yo = 
r;i = r;2 = to be the correct center and shape for this 
galaxy if the likelihood £ of the fit of p3l) to the observed 
data satisfies 

1^ = (14) 

and likewise for derivatives with respect to yo,rj\, and 772 . 

However we can also view the five functions in (|13p 
as the first elements of some complete basis set, with the 
coefficients of the last four terms serving as our centroid 
and roundness null tests. The basis-set method will consider 
the galaxy to be round and centered in this frame if the 
likelihood of the fit of (|13p to the data is maximized by 
nulling the coefficients of the last four terms. This is exactly 
the same condition for shape assignment as in the circular- 
model-fitting technique. Hence we can consider the latter 
technique to be a special case of fitting with a truncated 
basis set. 

MBLB report ~ 20% errors in measurement 
of deVaucouleurs-profile elliptical galaxies using Gauss- 
Laguerre measures of ellipticity — or worse with some other 
methods. MBLB report that these large errors may be at- 
tributable to a degenerate behavior of the fitting method: 
with very well-resolved deVaucouleurs-profile galaxies, the 
condition used to defined the size scale of the Gauss- 
Laguerre basis set may lead to attempts to measure shapes 
by fitting the central intensity spike of the galaxy. The basis 
functions can then be poorly resolved by the image pixel 
scale. We do not investigate in this paper the domain where 
the PSF-convolved image has structure well below the pixel 
scale, leaving this for future work. 



3 POTENTIAL SOLUTIONS 

The left panel of Figure [2] plots the bias in e measure- 
ments of Sersic-profile galaxies viewed at very high S/N ~ 
10* with a Gaussian PSF, as measured with the Gauss- 
Laguerre (GL) basis f u nction fitting method described in 
iNakaiima fc BernsteinI l|2007l . NB07). We see that galax- 
ies with low Sersic indices, close to the Gaussian value of 
n — 0.5, are measured with shape biases below the 1 part 
in 10'' threshold, as long as they are well resolved. However 
biases of up to several percent appear when the Sersic in- 
dex rises toward the deVaucouleurs value of n = 4, when 
the galaxies become highly elliptical, and when the galaxies 
are more poorly resolved. All three of these factors lead the 



galaxy to be more poorly represented by the finite GL basis 
set, opening the door to underfitting biases. 

One way to avoid the underfitting bias would be to find 
a basis set in which all galaxies are compact and represented 
wit hout signific a nt bia s by a finite number of coefficients, 
e.g. iNgan et all l|2009t) propose a basis set generated from 
Sersic profiles, which may represent Hubble-sequence galax- 
ies more compactly than the G auss-Laguerre basis sets used 
by NB07 and .Refregier fc Bacon ( 2003). A viable methodol- 
ogy from a Sersic basis has not yet been reported, however. 
The less regular morphology of fainter galaxies may make 
it difficult to obtain a sufficiently general basis. We will not 
pursue this approach in this paper. 

A second strategy is to simply extend the basis set to 
higher order until the biases are reduced to desired level. The 
dashed lines in the left panel of Figure [2] show the impact 
of increasing the GL expansion from order N = 8 (45 basis 
functions) to A'^ = 12 (91 basis functions). The biases are 
indeed substantially reduced in most cases, but remain un- 
acceptably high for poorly resolved objects and for objects 
with higher n and high ellipticity. We note however that 
a price has been paid, particularly for the poorly resolved 
objects, in that the resultant uncertainties in the ellipticity 
have risen substantially (see lower row of plots). While this 
is not a problem for an object with S/N = 10*, it will be 
a problem for real observations, where S/N is expensive to 
obtain. 

Why does going to higher order lead to higher noise 
primarily in the poorer-resolution cases? Recall that in the 
basis-set-fitting method, four basis functions are designated 
as the centroid and roundness tests. Considering just one 
of these coefficients bi and its associated function tpi, we 
know that as the basis set becomes more complete, the least- 
squares solution will converge (at high S/N) to the desired 
limit 



d"fcV.4k)J,(k) 



(15) 



We express the integral in Fourier space because it is easier 
to understand the noise properties of convolved images in 
Fourier space. We also work in a coordinate system where 
the object is nearly round, so that we are testing for 77 = 
and we can for the moment omit the S,, matrices for clar- 
ity. For background-limited images, the observed image /o(k) 
has a uniform (white) noise spectrum N. However to retrieve 
the pre-seeing image la, we must divide by the modulation 
transfer function (MTF), which is the Fourier conjugate of 
the PSF: L(k) = /o(k)/f'(k). This means that the effective 
noise power on the pre-seeing image is A'"/|T'^(k)|, rising to 
infinity at higher k. We thus expect that as the basis be- 
comes complete, the measurement of the roundness test bi 
will acquire noise level 



Var(60 =N (Tk U,(k)/T(k) 



(16) 



If we have chosen a roundness test ipi(k) which goes to zero 
before the inverse MTF 1/T'(k) blows up, then the noise on 
bi will remain bounded as we increase the fit to high order. 
We can reduce bias by fitting to higher order without paying 
a stiff noise penalty only if we have chosen our roundness 
test to roll off more quickly than the MTF in fc-space. 

If, however, we are using a roundness test that extends 
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Figure 2. Biases and noise levels for three shape measurement methodologies. In each panel the upper row gives the fractional error 
|(5e|/e on the measured ellipticity of a Sersic-profile galaxy with input ellipticity e labelled on the horizontal axis. The curves are red, 
green, blue, and magenta for Sersic index n = 0.5, 1, 2, and 4. The left three panels use the Gauss-Laguerre basis-set fitting method of 
NB07, deconvolving the apparent galaxy ellipticity from a Gaussian PSF with ellipticity 0.14 and half-light radius of half, equal, or double 
that of the galaxy, as marked. The deconvolved galaxy shape is modelled with order N = 8 (solid) or Af = 12 (dashed) decomposition. 
The middle set uses the FDNT shape measurement method introduced in this paper, executing the Fourier transform via a "native" GL 
decomposition of order 12 (solid) or = 16 (dashed). The right-hand set of panels shows biases using the FDNT method with FT of 
the real-space data. The lower panel shows the noise level in the inferred ellipticity for each case, relative to ellipticity error for a circular 
Gaussian galaxy with the same half-light radius and background noise as each test case, and therefore indicate how each measurement is 
being degraded by the ellipticity and non-Gaussian galaxies. We see that the FDNT methods dramatically reduce the shape measurement 
bias and with noise levels that are lower and less affected by order of the expansion. Note many FDNT measurements have errors below 
the scale of the plot, of the order 10~^, at the limit of precision of these tests. 



to regions where 1/T is large, then we will find the error on 
bi climbing higher and higher as we use more complete basis 
sets. As illustrated in Figure [3l this is precisely what hap- 
pens when using a model-fitting technique on an unresolved 
galaxy. One would like to choose a basis function t/ii which 
is of similar size to the (pre-seeing) galaxy la so that the 
representation of this galaxy is compact in the basis. How- 
ever if the galaxy is smaller than the PSF, then in Fourier 
domain, the roundness test ipi will extend to higher k than 
the MTF T, and the roundness test noise will diverge as the 
basis set becomes complete. 

When we use a roundness test that extends to k values 
where the 1/T deconvolution term becomes large, we face a 
tradeoff between noise and bias: as we truncate the basis set 
to avoid fitting degeneracies at high k, we become subject to 
bias if the truncated basis does not properly extrapolate be- 
havior to high k. We conclude that underfitting bias stems, 
at root, from an attempt to measure galaxy shapes using 
information at high k which is not present in the data. 

The MTF is not the only reason why our knowledge 
of a galaxy's la is incomplete. Finite sampling can cause 
aliasing ambiguities. In real space, missing pixels, and the 
necessarily finite measurement aperture for each galaxy, are 
also forms of missing information. We rely upon a finite-basis 
representation of the galaxy to fill in the missing information 
during a shape measurement. In the absence of a known 
finite basis representation of the galaxy population, our best 
strategy is to select a roundness test that is least sensitive 
to poorly observed regions of a:: or fc space. 



4 FOURIER-DOMAIN NULL TEST METHOD 
4.1 Designing the best roundness tests 

We propose a new shape-measurement method, "Fourier- 
domain null testing" (FDNT), which explicitly chooses a 
roundness test that does not use information destroyed by 
the MTF. We work in the Fourier domain where the effect 
of the PSF is more transparent. We assume that the round- 
ness test is linear and a pure quadrupole. Working in the 
coordinate frame where we are testing for 77 = the tests 
are: 

tlih) = j (fkia{\<i)w{k){kl^kl) (17) 

= J d^k [L{k)/f{k)]w{k){kl-kl) (18) 

t2{Ia) = j <fkia{\i)w{k)i2k^ky). 

Equation {TS} can be used to estimate the variance of ti 
in the case of white noise in the image. The response of 
the roundness test to a small shear is our desired signal; 
BJ02 show that when the galaxy is nearly round, this will 
be related to the unlensed intensity 7„ by 

|lcx Jd\kMk)k^{Ul.))^ (19) 

The brackets indicate averaging around the azimuth an- 
gle. Combining the expressions for signal and noise of ti we 
can find the radial weight 'w{k) that produces the Fourier- 
domain quadrupole roundness test with the best S/N on 
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Figure 3. Illustration of the difficulty of model fitting to a poorly 
resolved galaxy. If a galaxy is smaller than the PSF, then as shown 
here its Fourier transform extends past the MTF, which is the 
Fourier transform of the PSF. If we try to measure the pre-seeing 
shape of the image, we must divide the observed image by the 
MTF, which raises the power spectrum by the inverse square of 
the MTF, noted by the upper line. In the gray region of k space 
the image is very noisy. If we are testing galaxy roundness with 
a filter that is matched to the galaxy size, we will obtain a very 
noisy answer because the filter extends into the degraded zone. 



shear: 



w{k) = w,{k)To'{k), 



T(k)| 
d 



k^w,{k) = k—{L{k)) 



(20) 
(21) 

(22) 



The Wg term in the weight function serves to weight those 
annuU in k space where we expect the galaxy to exhibit 
the largest signal under shear. For a Gaussian galaxy with 



/a oc e 



the roundness-test weight should also be a 
BJ02 show that this choice 



Gaussian, Wg(fc) ; 
is also near optimal for an exponential-disk galaxy (Sersic 
n = 1). We will retain the Gaussian form of Wg even though 
it is suboptimal for steeper-cored Sersic profiles, because the 
details are not critical and any well-behaved function yields 
a valid roundness test. 

More critical is the T'o(fc) term in the weight, which 
serves to roll off the weight wherever the MTF is getting 
small. Note that the roundness test goes to zero at any ra- 
dius where the MTF drops to zero, so our roundness test 
never requires extrapolation into unmeasured regimes. As- 
suming that we have access to a Fourier representation of 
our observed galaxy, we do not need to fit the galaxy to any 



basis set to execute the roundness test. The tests for cen- 
troiding and roundness at shear rj are 



t. 



[d;,(lS„k|)/f(k)] J„(k) 



Wj{k) 



Wg{k)fi{k) X 



k^ - 


j 


= 1 


'2hx ky 


j 


= 2 


ikx 


j 


= X 




j 


= y 



(23) 



To assign an ellipticity to the galaxy, we first transform the 
observed galaxy /o(x) into Fourier domain, as well as the 
local PSF. For a given trial shear (771,772) we can readily 
create the four bracketed functions in (|23p . to evaluate the 
four null tests. Then we iterate the coordinate origin and 
the shear components (771,772) until all the tests are zero. 
We assign the galaxy an ellipticity e = tanh 77 at the null. 

The right-hand panel of Figure |2] shows the results of 
applying FDNT shape measurements to the same perfectly 
elliptical Sersic-profile galaxies measured with the GL basis- 
set technique. There is dramatic improvement: in nearly all 
cases the biases are too low to even appear on our plots. 



4.2 MTF shrink factor 

Note in the lower panels of Figure [2] that more elliptical 
galaxies have higher shape noise levels (expressed as errors 
in 77) , i. e. are less sensitive to applied shear. There is a subtle 
reason for this. The FDNT weight function contains a cir- 



cularized MTF term To{k) defined in Equation This 
averaging must be done in the sheared coordinate system 
where we are testing for roundness. Even if the PSF is cir- 
cular on the observed image, it will not be circular in the 
sheared coordinates. Furthermore the function Tq as defined 
in Equation (|2ip will change subtly as we test for roundness 
at different 77 values. In particular, if the MTF goes to zero 
at some fcmax in the observed frame, then the function Tq 
should go to zero at a smaller radius fcmax/s, where s is a 
"shrink factor" given by s = exp(77/2) = [(1 -I- e)/(l — e)]^''''. 
As illustrated in Figure |4j shearing the coordinate system 
makes the PSF broader in one dimension, effectively degrad- 
ing resolution in this direction. But since our roundness test 
requires a circularly symmetric weight function, we must 
shrink the entire weight function to fit it into the useful re- 
gion of the now-elliptical MTF. Because there is now less k 
space to use for the roundness tests, the shape measurement 
becomes noisier. 

Rather than recompute the circularized Tq via (|21|l at 
every shear being tested, we simplify by computing and tab- 
ulating a function T(fc) defined to be Tq in the observed 
coordinate frame (77 = 0). Then for non-zero 77 we just cal- 
culate the shrink factor s and use the function T(sfc) for the 
PSF part of the weight. 

A second effect of the distortion of the MTF as we 
change 77 is that the roundness tests are not invariant under 
shear of the underlying galaxy. A galaxy that is sheared to 
become more elliptical will be nulled by a roundness test that 
is slightly smaller (in k space) than the unsheared galaxy. 
This violates one of our original criteria for roundness tests 
and leads to potential multiplicative shear biases when the 
source galaxies have radial ellipticity gradients, a problem 
which we elaborate and solve in !j5] These problems are 
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Figure 4. Illustration of PSF "shrinkage" in k space for elliptical 
galaxies. We schematically draw the size of a galaxy and PSF in 
real space, in the Fourier domain, and in Fourier domain after 
shearing the coordinates to make the galaxy appear round. It is in 
this last frame that we want to apply a circular "roundness test" 
weight function in order to have well-behaved shape assignments. 
In the bottom row we see how a more elliptical galaxy forces the 
shaded roundness test to be smaller (in k space) if it is to remain 
inside the MTF where the data are not degraded. 



not specific to the FDNT method; all practical shape mea- 
surements apply some kind of PSF-dependent weighting to 
galaxies, and this weight will vary under shear unless we 
specifically design the weighting scheme to avoid this. 



4.3 Implementation notes 

The FDNT measurement can be quite fast, requiring an 
FFT of the pixelized la into Fourier domain, and the evalu- 
ation of the four weight functions at each point in the trans- 
formed image. The derivatives of the null tests with respect 
to centroid position are easily calculated; we solve for the 
centroid position by assuming linear response to centroid 
shifts relative to the starting position. The search for the 
shear that then nulls ti and t2 can use standard nonlin- 
ear root -finding techniques, in our case P owell's "dogleg" 
method (|Madsen. Nielsen, fc Tinglefj|2004 ). 

If the observed image has stationary noise, then the co- 
variance matrix of the Fourier components is diagonal, and 
it is also very fast to propagate the image noise into a covari- 
ance matrix for the 4 null tests. It is however computation- 
ally expensive to propagate a non-uniform noise map into 
Fourier domain, and therefore the FT method is ill-suited 
to estimating uncertainties on the null tests in this case. 



4.4 Fourier transform via basis sets 

Other shape-measurement methods make u se of Fourier- 
domain representations of galaxy images (|Zhanel l2008l : 
iMiller et aZ.ll2007l : iHosseini fc Bethg3 l2009l ). It may prove 
difficult, however, to turn such methods into practical weak- 
lensing algorithms, because the Fourier transform is not ro- 
bust to missing pixels, e.g. from cosmic rays. Furthermore, 
galaxies have neighbors, so we have a limited region over 
which to perform a transform without running into neigh- 
bor galaxies, which places a limit on the fc-space resolution 
one can obtain. Zero-padding a small postage stamp can in- 
crease the resolution of the FT, but one is then assuming 



that the galaxy fiux drops identically to zero outside the 
postage stamp, which may induce a bias in the shape mea- 
surement. 

We need some way to interpolate in real space to miss- 
ing pixels or extrapolate to large radii. We can call upon 
fitting methods again to do this. If we fit the observed im- 
age to a finite basis set such that /o(x) = ^ bi^pi{x.), then 
in Fourier domain the galaxy is /o(k) = ^6i?/'i(k). The 
roundness (or centroid) test tj is 



Cov(tj, tfe) 



i\ [w,(|s^k|)/r(k)] Mk), 

l3jlCov{bl,bm)l3mk. 



(24) 



(25) 
(26) 



Thus we need fit the observed image /o(x) just once to ob- 
tain the bi and their covariance matrix, i.e. we project the 
pixel data onto a description bi and execute the tests on 
the bi. If bi is of lower dimension than the pixel count, we 
have made the propagation of errors easier than in the FT 
method, and we become robust to missing pixels. 

We have implemented a version of the FDNT tests 
which uses the elliptical GL functions as a basis set. These 
basis functions have the advantage of being their own 
Fourier transforms (up to a power of i) , which makes the GL 
basis fit a fast way of bringing the galaxy into the Fourier 
domain. The fit is robust to missing pixels if the order of 
the GL decomposition is not too high, and the GL functions 
are compact in real (and Fourier) space so they do not need 
very large postage stamps to fit successfully. 

The central panel of Figure [2] shows the result of exe- 
cuting the FDNT method with Fourier transforms done by 
the GL decomposition. We see that for Sersic indices n — 2 
and n = 4, the GL transform induces biases that are worse 
than our FDNT-FT method and exceed the desired 10"^ 
threshold when the source ellipticity gets high. These galax- 
ies are of course the ones that are least compactly repre- 
sented by the elliptical GL basis, and hence the execution 
of the Fourier transform via GL decomposition leads to bi- 
ases. We would find larger biases if the PSF were an Airy 
function instead of a Gaussian, because it is known that the 
GL decomposition requires very high order to reproduce the 
fc = cusp of the Airy MTF (NB07). We have once again 
the situation that the use of a finite basis set to fill in miss- 
ing galaxy information can lead to biases if the basis set is 
not well matched to the galaxies. 

However we should note that the performance of the 
FDNT-GL method, while worse than FDNT-FT, is much 
better than the performance of the NB07 implementation of 
the basis-set fitting method. Note that the biases are much 
smaller than for the NB07 method, and also the shear noise 
level remains stable when we increase the order of the GL 
decomposition from 12 to 16, even for poorly resolved galax- 
ies, and the noise levels are generally lower. The FDNT-GL 
method differs from basis-set-fitting in that we are fitting 
the observed object profile, not attempting to fit the pre- 
seeing image. This makes the fit more stable so we can go 
to higher order. We can also use a basis that is matched in 
size and shape to the object, making the galaxy more com- 
pact in the basis. And by using the FDNT roundness test 
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instead of one of the basis coefficients, we make our shape 
measurement immune to noise that comes from attempting 
to fit the pre-seeing galaxy at k values suppressed by the 
MTF. 

We are not constrained to using the GL basis set; one 
can execute the Fourier transform by fitting any desired ba- 
sis set to the observed image. The direct FT method has 
the advantage of lower biases (and faster execution) than 
the GL FDNT method. However the benefit of projecting 
the pixel data onto a basis set before transforming is that 
we can compress the pixel data to a smaller number of pa- 
rameters, making propagation of errors more tractable, and 
also making the measurement robust to missing pixels or 
smaller measurement apertures. 



5 ELLIPTICITY GRADIENTS 

VB find that a galaxy composed of a superposition of bulge 
and disk components with different ellipticities was mea- 
sured with 2% multiplicative shear error, and that this er- 
ror was not reduced by increasing the radial freedom of their 
circular-model-fitting method (see their Figure 7). This mul- 
tiplicative error is a consequence of two circumstances. First, 
the galaxy's measured ellipticity varies with the scale of the 
weight function applied to the calculation of its quadratic 
moments, i.e. its core is rounder than its outskirts, or vice- 
versa in Fourier space. Second, the effective weight function 
of the circular-model-fitting method changes implicitly with 
the ellipticity 77 being applied to the model as the MTF gets 
sheared in the frame where the model is circular. As dis- 
cussed above, this violates the prescription that our round- 
ness test should be invariant under applied shear. We reiter- 
ate that this implicit change in the shape weighting function 
under shear is present in nearly every shape-measurement 
methodology published to date, with the choice of weight 
scale often buried deep in the implementation details. The 
resultant biases are probably endemic to current shape mea- 
surements. Here we investigate the problem and its solution 
in more detail within the context of the FDNT method. 



5.1 Estimation of bias 

Consider an unlensed, pre-seeing galaxy that when measured 
with a null test of scale s yields an ellipticity {rji ^ 0, r;2 =0). 
If the galaxy is not perfectly elliptical, it may have nonzero 
derivatives dr]\/da and drj2/ds with respect to a change in 
the weight scale factor. We will say that such a galaxy has 
"ellipticity gradients" even though its isophotes need not be 
ellipses. 

Assume that the weight scale s applied during a null test 
depends upon the amplitude T^t of the shear at which we're 
looking for the null. In most function-fitting methods, this 
change in scale size is implicit and not controlled, whereas 
we will control this explicitly in the FDNT method. There is 
some potentially non-zero derivative ds/drjt in most shape- 
measurement methods. 

Now assume that the shrink factor was s(r;i) when the 
galaxy was unlensed, and then we apply a small lensing shear 
(Aryi, A772) to the galaxy and we remeasure it. The ellipticity 



(77^,772) measured by the method will satisfy 



m 

As 



771 + ) + A771 



ds ) ^ tanhj7i 
ds 



Ar;2 



(27) 
(28) 
(29) 



The key observation is that the parenthesized quantities in 
Equations (|27|l and (|28p are the shape of the unlensed galaxy 
if measured with the slightly different scale s + As that will 
be used when measuring the lensed galaxy. The right-most 
terms are the differential to the ellipticity caused by the ap- 
plied lensing shear, as per Equation Qlj The applied shear 
shifts the measured shear according to 



5ri2 



Aril 



1 



dr}i ds 
ds drjt 



-1 



■A772 -I- -3— 07?1. 

ds drjt 



tanh 771 



(30) 
(31) 



The ellipticity gradient hence causes the measured response 
to shear to depart slightly from the applied shear. Be- 
cause the coordinate system was chosen to align with the 
galaxy (772 = 0), a change in 772 is a rotation of the 
ellipticity — we would find the change in position angle 
dp/ds = (2/771)^772/^5. This is a quantity that changes sign 
under reflection of the unlensed galaxy about the x axis. 
Since the unlensed population of galaxies should be parity- 
invariant, we can assume that {dri2/ds) — when averaged 
over the population, so the ellipticity-gradient effect on ^772 
averages to zero. 

The ellipticity gradient therefore causes a rescaling of 
the measured shear 577|| component lying parallel to the 
galaxy's major axis, while leaving the orthogonal shear com- 
ponent Srj± unchanged on average. To gauge the effect on 
the weak lensing shear measurement, we average the mea- 
sured differential shear over a population of galaxies with 
unlensed ellipticity 77 and random position angles /3. After 
some algebra to execute the average over /3 we obtain, for 
both shear components j £ {1, 2}, 



7^l 

EG 



= TZiArj-j 

_ 1 

^ 2 



1-EG 
d77|| ds 
ds drjt 



+ 



tanh 77 



(32) 
(33) 

(34) 



5.2 Correction for ellipticity gradients in FDNT 

In the FDNT method, the weight factor is constructed 
explicitly from two fc-domain factors, the galaxy function 
Wg{k) and the circularized PSF factor T^{k). If either 
changes when the galaxy is sheared, we will incur an elliptic- 
ity gradient bias. We have already described that we replace 
the Tq term with T(fcs) with s — exp{r]t/2) to account for 
shrinkage of the MTF minor axis in the circularizing coor- 
dinate frame. We therefore have ds/drjt ~ exp(77t/2)/2. 



^ Note that we are describing the galaxy ellipticities with rj' in- 
stead of the e' = tanh??' used in Equation (|4j. 



10 Bernstein 



We also need to control any change in Wg (k) under shear 
of the galaxy. We first set Wg{k) = Wg{kag) where Wg is a 
fixed function form — a Gaussian — with explicit scale ag . We 
choose ag by requiring solution of this equation: 



(f'kk'^T{ks)Wg(kag) 



(fkk^Tiks) [Jo(S,7^k)/f(S,7^k)] . 



(35) 



This condition requires that, in the sheared coordinate sys- 
tem which make the galaxy round, the fe^ moment of the 
weight function must match that of the galaxy itself. Both 
functions are weighted by the factor T(fcs) in order to con- 
fine the calculation to the domain where the measurement 
survives convolution with the PSF. Because the condition is 
applied in the circularizing coordinate frame, the chosen Gg 
should be invariant under applied shear at fixed s. 

The roundness-test weight is therefore scaled by s both 
directly through T{ks) and indirectly through any potential 
variation dog/ds implicit in Equation (|35p . We can calculate 
the resultant dr]\\/ds either by propagating derivatives of the 
roundness tests with respect to s and (jg , or by intentionally 
perturbing s and repeating the solution for the rj which nulls 
the tests. In particular, if t is the vector of null tests, then 
the measured ellipticity e of the galaxy varies with s as 



dr] 

ds 



dt 

ds 



(36) 



We can then take the component of this derivative in the 
direction parallel to the galaxy ellipticity to get dr]\\/da and 
derive the ellipticity-gradient response factor EG. 



5.3 Testing 

For galaxies that are not perfectly elliptical, there is no sin- 
gle correct measure of the galaxy ellipticity, and we must 
perform a "ring test" as described by NB07 in order to test 
for shear measurement biases. Copies of the unlensed galaxy 
are drawn with rotation 9 uniformly distributed from 0-27r. 
Each copy is then sheared by an amount Ar;, taken here 
to be along the x axis. The copies are then convolved with 
the PSF, sampled on the pixel grid, and given appropriate 
noise levels. Each copy is then measured using the FDNT- 
FT method to yield galaxy shape components 771, i and 772,1 
and an estimate of the ellipticity-gradient factor EG. From 
EG and the measured shape 77 of each galaxy we also obtain 
an estimate TZi^i of the shear responsivity in Equation (|33|l . 
From the ring ensemble we then estimate the applied shear 
as 



A?7 : 



(37) 



If this estimator reproduces the shear A?7 applied to the ring 
without bias, then it will also produce an unbiased shear 
estimate from any source population that is isotropic before 
lensing. When the ring test is executed with applied shear 
Ar7 = 0, we always find the estimator consistent with zero, 
indicating that we have no additive biases arising from the 
non-circular PSF. We search for multiplicative bias by seeing 
how well an applied shear of A?7 = 0.02 is retrieved. 



We do not report the mean of the estimated component 
772 orthogonal to the applied shear. Symmetry considerations 
imply that the mean 772 will vanish when the source popula- 
tion is invariant under reflection about the x axis — and this 
does occur with ring test for Test Galaxies I and 2. Test 
Galaxy 3 does not have the reflection symmetry, so elliptic- 
ity gradients or other biases could cause (772) 7^ 0. However 
any Universe containing galaxies like this would contain an 
equal number of mirror-imaged versions of Test Galaxy 3. 
As argued in i]5.1l this will cancel 772 biases from ellipticity 
gradients. 

We note that the "perfect" ellipticity transformation 
Equation Q and the resultant responsivity estimate (|33|l 
take a linearized response of the galaxy shape to an applied 
shear. For any measure of shear, there are also terms cubic 
in the shear which are significant if we have applied shear 
A77 ~ 0.05 and wish to measure Ar; to < I part in 10^. In the 
Appendix we give the third-order responsivity TZz that must 
be included to measure shears below the part-per-thousand 
level. The results presented in this section apply the third- 
order corrections expected for a perfect shape measurement 
method to the ring test results. The cubic terms derived in 
the Appendix also include transformation of the estimated 
shear from the unit-determinant A?7 to the traceless 7 mea- 
sure. 

We first perform a ring test using the FDNT-FT method 
on the simulated spiral galaxy of VB. This galaxy has 20% 
of its flux in a "bulge" with Sersic index n = 1.5 and ellip- 
ticity e = 0.05 and half-light radius of 1.5 pixels. The "disk" 
component has n = 1, e = 0.2, and half-light radius 11.8 pix- 
els, with the same position angle as the bulge component. 
This galaxy is drawn in the top left panel of Figure [5] and its 
measured ellipticity is plotted in the panel below. The mea- 
sured ellipticity changes as we alter the size of the applied 
PSF and hence the scale of the null test weight function. 
The dashed line in the lower panel shows the multiplicative 
error on the applied shear as measured by the ring test with 
no correction for the ellipticity gradient. The star marks the 
2% bias of the VB method, slightly more biased than our 
uncorrected FDNT method. If, however, we calculate the el- 
lipticity gradient correction EG (third panel) and apply it to 
the responsivity as in Equation H33p . the error in the inferred 
shear is reduced 100-fold to negligible levels of < 3 x 10"'*, 
as recorded in the bottom solid line. 

We conduct tougher tests by producing galaxies with 
larger ellipticity gradient corrections. In the middle panel of 
Figure[5]the two Sersic components have orthogonal elliptic- 
ities and the required EG correction is up to 4 x larger than 
for the VB galaxy, yet the i5G-corrected shear estimate is 
still accurate to < 3 x 10~^ . Test Galaxy 3 in the right panel 
is a true torture test: the two components are not orthogo- 
nal and not concentric, so the galaxy has multiple maxima 
has no inversion or reflection symmetries, and has isophotal 
twists. Without EG correction, the multiplicative errors on 
shear are as large as 6%, but the calculated EG reduces bias 
at least 20-fold, close to the 10""^ target. 

We speculate that the residual multiplicative errors 
in the torture-test are due to the limitations of linearized 
analysis of the ellipticity gradient effects in deriving Equa- 
tion dSni. For the VB galaxy, EG « 0.03, while EG ~ 0.1 
for the torture-test galaxy. We should not be surprised to 
find shear errors at 0{EG^-y) or 0{EGj^). The the resid- 
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Figure 5. Multiplicative errors in the shape measurement of 
a simulated galaxies, having an ellipticity gradient, using the 
"FDNT-FT" method described in the text. Each test galaxy 
has an outer exponential disk and an inner deVaucouleurs-profile 
"bulge," with the two components having different ellipticities 
and scale lengths so that the overall ellipticity of the galaxy de- 
pends on the scale of the measurement weight. For each test 
galaxy, the second row plots the ellipticity measured vs half-light 
radius of the PSF. The third row plots the factor EG, measured 
from the data, which should be applied to the FDNT-FT shapes 
to infer the applied shear as per Equation H33| l. The bottom 
row plots the multiplicative error on the inferred lensing shear 
both with and without the applicatio n of the EG correction . Test 
Galaxy 1 is the "spiral" galaxy from lVoigt fc Bridlj | |2009|) . with 
a mild change in axis ratio, and the multiplicative errors are neg- 
ligible after EG correction. Test Galaxy 2 is a harsher test, with 
more elliptical and orthogonal components, and still has negligible 
error after correction. Test Galaxy 3 is a yet more difficult case, 
with multiple maxima, a larger ellipticity gradient, and no inver- 
sion symmetry — even in this case, the EG correction improves 
the multiplicative error 20-fold, to ~ 10~^. 



ual errors appear to indeed scale with applied shear as 7 , 
but further investigation is required to check dependence on 
EG and other galaxy or PSF properties. 

We note that "stacking" methods will not suffer from 
ellipticity gradients, since stacking creates a m ean galaxy im- 
age that ten ds toward circular at all isophotes l|Kuiikenll999l : 
lLewis|[2009l ). A stack of galaxies with random (unlensed) 
orientations has ellipticity \/iV smaller than the typical in- 
dividual galaxy. Hence the ellipticity gradient of the stack 
is reduced by y/N as well, as is the induced EG bias. For 
N > 1000, stacking will reduce the ellipticity-gradient bias 
as well as the EG correction presented herein. The relative 
success of stacking methods in the GREAT08 tests is likely 



a reflection of this advantage over other previous methods. 
The FDNT method is the first proposal to properly correct 
for the ellipticity-gradient bias without stacking. 



5.4 GREAT08 tests 

We apply the FDNT-FT measurement method to test im- 
ages comprising the low-noise sectio ns of the GREAT08 
Challenge. As extensively described in lBridle et al.l (|2009ah 
and Bridle et al. (2009b), this challenge was to estimate the 
lensing shear applied to each of 15 images, with 10** sim- 
ulated galaxies of S/N ^ 200 in each image|f] GREAT08 
results are scored by Q = lO^^/a^^MSi where cr ^.^rms is the 
mean error in the estimate of each component of shear. 

We found that the postage stamps on which the 
GREAT08 galaxies were drawn were small enough to trun- 
cate many of the more elliptical test galaxies, and also too 
small to contain the real-space versions of the filters being 
selected by the FDNT-FT algorithms. D. Kirk and S. Balan 
kindly produced a new realization of the low-noise blind 
GREAT08 test suite, drawn from the same galaxy distribu- 
tions as the original, but placed onto larger postage stamps. 
The FDNT-FT algorithm was not provided with any knowl- 
edge of the characteristics of the underlying galaxy popula- 
tion, and we were blind to the input shears for the "large- 
stamp" test images. We did, however, take advantage of the 
exact analytic form provided for the PSF. 

The be s t perf ormance repor ted for any method in 
iBridle et af] (|2009b') was Q ^ 500. iGruen et all (|2010l ') have 
subsequently reported Q ~ 1500 for some subsets of the 
GREAT08 images, using a neural network to calibrate shear 
measurements to simulations of known input shear. Without 
any training or recalibration, the FDNT-FT measurements 
achieved Q = 2997, corresponding to an RMS shear esti- 
mation error of 0.0018, which is only 20% higher than the 
expected statistical uncertainty due to noise in the input 
ima ges. This Q value i s 6x higher than any other reported 
m Bridle et alj (|2009bl ). and 3x higher than the value es- 
timated for the shear accuracy required to have negligible 
degradation of the statistical errors of hemisphere-scale WL 
surveys. 



6 CONCLUSIONS 

Seeking a weak lensing shear measurement methodology 
that can achieve systematic errors below 1 part in 10^, we 
examine the origin of the larger biases reported by VB and 
MBLB for measurement of noiseless galaxies using several 
current object-fitting methodologies. We show that shape 
biases can be generated when fitting with the basis sets are 
(necessarily) truncated to finite size or with perfectly ellip- 
tical galaxy models. A multiplicative shear bias will typi- 
cally occur when the radial profile of the galaxy is not well 
matched by the finite basis set. This could be remedied 
by choosing a finite basis that is complete over all galaxy 
shapes, but such a basis has not yet been devised. 

One could extend the fitting functions to higher order to 



^ See'http : //www . greatchallenges . info for documentation and 
images from GREAT08 and its successors. 
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reduce the bias to acceptable levels — VB show this to work 
for most of their cases — however we demonstrate that there 
is a fundamental limitation to this approach. If the "round- 
ness test" being applied to the pre-seeing image extends to 
regions of k space that have been highly or fully suppressed 
by the MTF of the image, then any attempt to fully elimi- 
nate the bias will cause the noise in the shape determination 
to diverge. This will occur whenever one is trying to use a 
basis set that is size-matched to a pre-seeing galaxy image 
that is smaller than the PSF. 

We propose a new shape measurement, "Fourier- 
Domain Null Tests" (FDNT), which constructs linear func- 
tions of the galaxy's Fourier representation that are explic- 
itly designed to maintain bounded noise in the presence of 
an MTF, and which are pure dipole and quadrupole func- 
tions of the pre-seeing galaxy intensity pattern. Following 
the BJ02 and NB07 prescriptions for "geometric" shape 
measurements, the centroid and shear of the FDNT func- 
tions are varied until the tests are nulled, and the galaxy 
is assigned a shape according to the shear that nulled the 
tests. This shape is guaranteed to transform in a well-known 
way under applied lensing shear. 

No model fitting is required for the FDNT method if a 
full real-space observed image is available. If the image has 
missing pixels or a restricted mask, then we can use basis- 
set-fitting methods on the observed image to interpolate the 
missing data and estimate the Fourier-domain image. We 
demonstrate that the FDNT method obtains negligible mul- 
tiplicative errors when applied to galaxies with pure ellipti- 
cal isophotes, even for poorly resolved deVaucouleurs-profile 
galaxies with high ellipticity. When using basis sets to es- 
timate the Fourier transform, measurable biases return, de- 
pending upon the choice of basis set, but these biases are al- 
ways smaller with FDNT tests than with the NB07 method. 

We identify a second shear-measurement bias that is 
likely even more pervasive. This one occurs when the mea- 
sured ellipticity of the underlying galaxy varies with the 
scale of weight function applied by the measurement. When 
we measure the shape of such an "ellipticity gradient" galaxy 
we must be sure that the weighting function stays fixed if 
the galaxy is sheared. In forward-fitting shape measurement 
methods, the effective weighting function arises implicitly 
from a combination of the fitted galax;y model, the noise dis- 
tribution, and the PSF. This implicit function does typically 
change as a galaxy is distorted, inducing a shear measure- 
ment bias. Since our FDNT method controls the weighting 
function explicitly, we can determine and correct the elliptic- 
ity gradient effect on shear measurement. We demonstrate 
that this reduces multiplicative biases by factors of 20-100 
for test galaxies with ellipticity gradients. The result is a 
shear measurement method yielding multiplicative biases in 
ring tests that in the worst case are at the target 10"'^ level, 
and well below in most cases. 

Both of the biases that we describe arise because trun- 
cated basis-set-fitting methods are weighting the pre-seeing 
images in a way that we do not explicitly control. The FDNT 
method makes the weighting of our roundness tests explicit 
so we can avoid degraded regions of k space and correct for 
ellipticity gradients. One might be chagrined that the EG 
corrections are an empirical calibration correction that must 
be measured on the (noisy) galaxy images, much like the 
polarizability factors of the commonly used KSB method 



(|Kaiser. Squires. &: Broadhurst I Il995l ). Certainly noise on 
the estimate of EG is an issue, since it enters non-linearly 
as 1/(1 — EG) in the calculation of responsivity. We note 
however that polarizability corrections are typically 0(0. 1) 
or larger while the EG correction is this large only for our 
"torture test" galaxies, and would more typically be 0.01- 
0.03. Since the EG correction are smaller to start with, they 
do not need to be determined to as high relative accuracy 
as the KSB polarizabilities. 

We have investigated an idealized case in this pa- 
per, with well-sampled images having no missing pixels, no 
crowding, uniform white noise, Gaussian PSFs, and very 
high S/N on a single image. Even in this case we note that 
obtaining < 10"'' multiplicative errors requires great care, 
e.g. the use of nonlinear shear responsivities as described in 
the Appendix, and knowledge of the PSF to very high ac- 
curacy. In a future paper we will extend the FDNT method 
to finite noise. We note that extension of the FDNT method 
to measurement of a galaxy imaged in multiple exposures is 
straightforward — one simply creates a null test that is the 
5'/A''-weighted sum of the null tests for the individual expo- 
sures. We iterate the common centroid and shape of these 
tests until the summed test is nulled. Each exposure can 
have its own PSF, and unlike fitting methods, we need not 
even assume that the galaxy has the same underlying ap- 
pearance in each exposure. So FDNT is easily applied to 
the full collection of multicolor images collected for a single 
galaxy. 

This work is supported by grant AST-0607667 from 
the NSF and DOE grant DE-FG02-95ER40893. We thank 
the anonymous referee for his/her help in clarifying the 
paper. We also thank Donnacha Kirk and Sree Balan for 
their substantial efforts in producing new realizations of 
the GREAT08 low-noise images for the blind tests reported 
herein. 
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APPENDIX A: THIRD-ORDER SHEAR 
RESPONSIVITY 

Consider a population of A'^ perfectly elliptical galaxies with 
axis ratios a/b — tanh(j7/2) and position angles uniformly 
distributed around < /3 < 27r. If these are distorted by a 
lensing field with shear 7 such that they are magnified by 
1+7 along the x axis and demagnified by 1 — 7 on the y 
axis, then a measurement of the lensed galaxies will yield 
ellipticities rj' and position angles /?'. We assume that we 
will infer the shear by calculating the mean of the ellipticity 
component r]i: 



^'^1^ ^J^^ cos 2/3'. 



(Al) 



Here we allow for application of some weighting function 
w{ri') to the averaging, as suggested in BJ02 to reduce shot 
noise in the estimator. We wish to determine how (77^) de- 
pends upon the input shear 7, keeping terms to third order 
in 7. The exact equations for shear addition are equations 
(2.12) of BJ02. In our case they become 

(1 + 27^) cosh rj + 27(1 + 7^) sinh rj cos 2/3 
sinhr;sin2/3 (A2) 



cosh 77' 
sinh 77' sin 2/3' 



After solving these equations for r]\ in terms of 7 (to third 
order), 77, and /?, and then averaging over /3, we obtain 



7^3(^) 



277^l(^7) + 27^7^3(77), 

■qw + TO + Yw, 

y2(y- 1) 



(A3) 
(A4) 



Y + 



+2r]w' 



77^ 
2Y{Y 



77^ 



+ 



3Y^ 
2^ 



Y 



+w"(2y + 3) +77TO" 



(A5) 

, • (A6) 
tanh 77 

When observing a population of galaxies, we will not know 
the unlensed 77 values, so the responsivities TZi and TZ^ will 
need to be calculated at the observed (lensed) values 77'. The 
relation of TZ\ as determined from the lensed distribution to 
the TZi at the unlensed 7; is, to second order in 7, 



^ 7^l (77') = iv [7^l (7?) + 27'7^3 (77)] 



(A7) 



Por the cubic term TZz we need only note that the mean 
value is the same, to linear order in 7, whether we average 
over lensed or unlensed 77 values. 

We assume now that we have a catalog of observed el- 
lipticities 771, i, for galaxies i G {1, . . . , A*'}, each assigned a 
weight Wi and the linear and cubic responsivities TZi^i and 
7?.3,i assigned from the observed ellipticity for each galaxy. 
We define the linearized shear estimator 



71 



2Y.n1,: 



(A8) 



then the above equations imply that the expectation value 
of the linearized shear estimator will be related to the true 
reduced shear via 



(71) = 71 1 - 7 



(A9) 



The same equation holds for the other shear component. 
Given an linearized estimated shear (71,72), one can quickly 
solve for the true shear (71,72) that would produce this lin- 
earized estimate. Alternatively the power spectra, etc., of 
the linearized estimator can be predicted from theoretical 
models for the power spectrum (and higher-order functions) 
of 7. 

We validate numerically that this procedure yields shear 
estimates accurate to 0(7^) for any shape measurements 
which transform under shear according to the exact shear 
addition formulae from BJ02. 

In this work we show that most shape measurement 
methods have an ellipticity gradient correction EG that al- 
ters the mean shear response along the principal axis. The 
expression for TZi should be altered from Equation (|A4|I as 



7^l(T7) = 



rjw + w 
1-EG 



+ Yw. 



(AlO) 



We have not calculated the effect of the ellipticity gradient 
on the cubic responsivity TJ-a. 



